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THE COMPUTATION OF 77 

PREFACE 

Althuugh the familiar »ymiH)l in) fur pi did not come into general use until a 
little over two hundred yearsi ago, computing the numerical value of tt has engaged 
the attention of mathematicians from the time of the ancient Egyptians down to 
the electronic computers of tiniay. Thus the Ahmcs Papyrus (Egypt) of about 
1800 BX. gives the area of a circle as 



HI 



where is the diameter. Thit$ is equivalent to taking tt as 



or approximately 3.1604 ... At about the same time, the Babyhmians, the Hindus 
and the Chinese took tt as equal to 

The early Greeks were c<mcerned with the problem of squaring the circle, and 
in the course of his searching, Archimedes, about 250 B.C., assumed the value of 
ir to lie between 

(=»3J408...) and3|(= 3,1428..,). 
The (Chinese scmn decided Ubout 100 A.D.) that tr was approximately equal to 

vTo, 

or 3.162 . . , About 150 A.D.« the renowned Greek astronomer, Ptolemy of Alexandria* 
Using the sexagesimal system of notation, stated that ir — 3^'30^, or, as we would 
write it today. 

60 {60)3 

which gives the approximation 3J416, or 3.141666 . , . The Hindu mathematician 
Aryabhata, about 500 A.D., gave two values of 7r, 

^ 177 ^ 62,832 

"^255 ^'"'^ 20,000' 

both of which give the value 3.1416, exactly. The latter fraction is presumably 
cakulaled from the perimeter of an inscribed p<ilygon of 384 sides. 

For the m thousand years or more mathematicians in many lands struggled 
with the problem^ but with little progress. Finally, about 1580, Francois Vleta, a 
pioneer French algebraist, using a polygon of 393,216 sides, found tr correc t to 

• • « 

111 
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nine decimal places, placing it between 3J4IS926535 and 3.14I5926S37. It appears 
thai be was the first mathematician to use an infinite priniuct in this connection, 
asserting that 

By way of contrast, contemporary mathematicians have computed the value of 
7r to more than 100,000 decimal places* This is a fantastic achievement from any 
point of view. The story of this long evolution lasting nearly 4000 yeai^ is indeed 
fascinating; the ciima% is vividly set forth in the present collection of essays. 

William L. Schaaf 
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FOREWORD 

The numerical value ir can be approximated by either of two general methods 
with ai close an approximation to its *Unie'* value aa we wtih. One method m 
geometrical. Thin is the classical approach first used by the Greek geometers and 
by mathematicians generally until comparatively modem times, that is until about 
16S0. It involves computing the perimeters of (Kilygons inscribed in and circum- 
scribed about a circle, and assuming that the circumference is intermediate 
between these perimeters. As the number of sides of the polygons is increased, 
the approximation becomes more accurate* In fact, if the areas of the polygons 
are used instead of the perimeters, an even letter approximation can be obtained. 

The second method, the modern approach, depends upon an expansion of tt into 
some equivalent analytical expression such as a converging infinite series or a 
convergent infinite product. One of the first mathematicians to use such an expres* 
ston was Vieta, as we have already seen. Another was John WalUs, who showed, 
in 16S6, that u could be expanded into the infinite pnnluct: 

n 2244668 
2*1 '3'3'5*5'7*7 

Many other nmthematirians have developed various expansions for evaluating 
IT, among them James Gregory. G. W. von Leibniz, John Machin, Leonard Euler, 
and C. F. Gauss. 

Perhaps a few words of explanation about infinite series will help you to under* 
stand the following articles better. 

A succession of numbers which follows a definite law or pattern is called a 
finite seqi^ence; for example, 

(a) 2, 4, 8, . . . 2% 

run ^ 1 J 
or (b) U^.g^.-.g^- 

where n is a ptjsitive integer. 

A sequence that is endless, having a first term but no last term, is called an 
i^finiie sequence* 

If we consider the sum of the first n terms of a finite sequence, we refer to the 
indicated sum as a finiie series; thus 

*i*{n*) ^ 1 +4 + 9+ 16 + 25 + 36=91. 
»- f 

We designate the indtrated sum of the terms of an infinite sequence as an 
ir\^nite series. But this is not a sum in the usual sense of the wc»rd, because the 
terms of an infinite series can never all be added term by term* 

1 
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If the nuci'eiiiihm of partml nuttu uf an infiiittc* Mcrii'tt inrreaMi^ii iitdefinifeiy ai) 
n iticreaiipp indefinitt'lv, ihi* st-fit-s in ^aid to bo diverfi^nt, and the •'sutir' uf tlu* 

If. nn the other hand, ihn MU(M*e^t«ioti of parlial mmn of an inKnitf Hi*rWh ap- 
liriNif lu'K a iimifing value as n invt^i^nvn tndrHnit«*l>, thr m-nvs is said tit bt* con* 
t^rgeni, and thf •'ikuni'* of tht* series refers to liiii^ iimitiug value. For example, 
in the inAnite series 

1+1+1+1+ +1+ 

the lintiting value is 2. The sum of any Htiile number of terms of this series, how- 
ever great, is always^ les:^ tiian 2: but by taking more and more terms, the partial 
sum ean be made as efusi* to 2 as ne wish. The limiting value "2" is ealled the sum 
of the eunvergent inhnite serit»s iti i|uesti<io. 



The Latest About n 



by HoHwd Ei>eM 

On July 29, 1961, Dr. Daniel Shankn and Dn John W. Wrench, Jr. cumputed 
IT to 100J265 decimal places un an IBM 7090 system in the IBM Data-center in 
New Yurk. The computation took 8 houn 43 minutes, including 42 minutes tu con- 
vert the final result from binary to decimal form. A check run, using a second 
formula, confirmed the accuracy of the first run to 70,695 decimals, and subsequent 
runs on 7090 computers in the Washington area showed that a machine error oc- 
curred in the initial run. Dr. Shanks and Dr. Wrench now have results that agree 
perfectly (including the convention and printing) to 333,075 bits or 100,265 decimal 
places. 

The first computation employed the formula 

TT - 24 iw Hi) + 8 tan U^) + 4 tan Mjss)^ 

which was published by Carl St%mer in 1896, This formula is especially well 
adapted to binary ciimputers, inasmuch as the evaluation of pi^wers of ^ on such 
computers can be accomplished simply by shifting. 
The check compulation was based on the formula 

48 tan-»(,^) + 32 lan U^) - 20 tan"'(^), 

of Gauss, which was used by George Felton to compute tt to 10,021 decimal places 
on a Pegasus in 1958. Because of the overlapping terms in the two formulas used, 
the check consisted almost entirely of the computation of 48 tan"*(^). This re* 
quired 4 hours 22 minutes on the IBM 7090 system. 

Following the discovery of a machine error in the evaluation of 24 tan~Hi)be- 
yond the 70«695 decimal place, all the arctangents were individually recomputed 
on a second 7090 system, and complete agreement was reached in all phases of 
the calculation to 100,265 places. 

On September 11, 1961, the 7090 system prepared a count of the frequency dis- 
tribution of the decimal di^ts of n in successive chiliads. Comparison with this 
lateftt count, carried to 100,000 places, revealed a few errors in Dr* Wrenches 
enumeration of the distribution of the digits 7, 8, 9 as published in Table 1 of his 
paper, "The evolution of extended decimal approximations to ^r/* in The Mathe* 
MATICS Teacher, LIII (Dec.. 1960), 648. This earUer count had been based on a 
computation of tr to 16,167 decimal plrres obtained on July 20, 1959, using a pro* 
gram of Francois (^nuys, on an IBM 704 system at the Commissariat a TEnergie 
Atomlque in Paris- In Table 1 of Dr. Wrench's article in The Mathematics 
TcachEH, the last four entries in the 7-coiumn should read 1258, 1342, 1439« and 



1546, respectively. In the S^olumn read 1243, 1336, 1455, and 1543, and in the 
9H!oIumn read 1306, 1418, 1513, 1615. With the»e correcUons, Table 1 is entirely 
free from erroni. 

Furthermore, on August 22, 1961, Dr. Shanks and Dr. Wrench also computed 
e to 100,265 decimal pUces on an IBM 7090 system in 2 hours 25 minutes, exclu- 
sive of the conversion to decimal form, which again required 42 minutes. The 
well-known factorial series was used, and a total of 25,266 reciprocal factorials 
were evaluated lo the stated accuracy. This confirms the 60,000 decimal place 
computation of e on the Illiac by D. J. Wheeler in I^^cemker, 1952. Wheeler*s 
calculation required 40 hours on the Illinois computer. The accuracy of the 100,265 
decimal places constituting thin latest approximation to e has been confirmed by 
a second calculation, which gave the respective sums of the even^- and odd* 
numbered terms of the factorial series, yielding approximations to both e and l/c 
to this accuracy. 

Dr. Shanks and Dr. Wrench have prepared a joint paper on their calculations of 
ir, which appears in the January, 1962, issue of Mcukematics of Compuiaiion. 
Appended lo their paper is the value of n truncated to 100,000 decimal places. 
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Forewofil 

As it turns out, tome infinite series converge more rapidly than others. 
Cooaider, for eEampIe, the aeries 

jc* 

arctanx=5s;i.^j + — ^jj 
If we set X 1, we gel Gregory's series, 

+ ; * (2) 

but tiiis series converges too slowly for purposea of computation. 
On the other hand. Machines formula, 

TT . 1 1 

- « 4 arc tan — arc tan — (3) 

used in conjunction with the expansion (1) alxive, converges much more rapidly. In 
fact, you can get a rather good approiumation to n simply by taking the first four 
terms of (1) when x ^ 1/5 (or .2), together with first term of (1) when x ^ 1/239. 
Try it and see for yourself! 

The reader who is familiar with trigonometry, may be interested in the deriva- 
tion of Machines formula* To prove that 



7 4 arc 
4 



Un(|)-«rctan(^), 



iet arc tan 1/5 = a, so that tan a 1/5. (1) 

2 

Ti. . o 2 tan o 5 5 

I hen tan 2a = = — = — = — . (2) 

1 - tan* a 1 JL 12 ^ ' 

25 

10 

and tan 4a = — — = — ^ (3) 
95^ 1 19 • ' ' 

144 

Since tan 4a is very nearly equal to 1, we see that 4a is approximately equal to ir/4. 
Now let 4a 7r/4 + arc tan x. (4) 



Recdl that tan + f ) - = |± 

V 4/ l-t&nA 1- 



tan A 



tan A 



Hence. ~ = tan 4a tan ^arc tan * + f ^ ~ • ' "° 2I9 

Therefore, from (4), we have 7 =" 4a — arc tan x, 

4 

or I - 4 arc tan (|) - arc tan . 



(5) 
(6) 
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A Series Useful in the Computation of tr 

by J, S. Frume 

One of the standards ways of cumputtng tt is based on Machines formula: 
and the series expansion 

(2) lttn-'i = *-^ + ^-... 

W. Shanks used precisely this in computing tt to 707 decimal places. In applying 
this aeries to the case x = 1/5, the individual terms are easily computed as deci- 
mals, and the series converges rapidly enough so that 55 terms suffice for 50-pIace 
accuracy. When we set x == 1/239, however, the individual terms, involving powers 
of 1/239, are not easily expressed as decimals, so that computation beyond 15 
decimals is laborious despite the rapid convergence. If, however, we expand the 
lerm^ in powers of 1/240, we obtain a new series which converges rapidly, and 
whose terms are easier to compute as decimals. The result is expressed by the 
formula: 

(3) 

The terms are alternately positive and negative in groups of three, so the error 
in breaking off the series is less in absolute value than the first group omitted. 
The series converges for |l| < 1/VX. Setting t = 1/240, we obtain the series 



I 1 2/1 \^ 2/1 V 
^^"^'239 ^ 246 ^2 1245/ 3 V246/ 

5 \240/ 6 V240/ 7 V240/ 



The computation is conveniently arranged as follows: Divide 1 by 240, this by 120, 
this in turn by 240, and so on alternately. This takes care of the numerators auto- 
matically. It remains only to divide each term by the corresponding exponent, and 
add and subtract appropriate terms. Sixteen terms give 50-place accuracy. 
The proof of formula (3) is a special case of the following: Let 

at^h ia c 

cl + a — ib + d 

7 
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Then 

ict + d) " Hat + b) 



ic + ia)t+id-^ib) . _d±ib^_^, d+ib 
"^'^ {c-iaU+id-'ib) ^ , , c-u», "^^"«rf^' 

^ d-ib' 

This series converges for |f { < 1/r, but it is useful for computation only when the 
values of sin a6 are convenient rational quantities. If 2 °« t, we have x^t, and 
obtain the series (2). The other case of interest is 2 = 1 + i, at * 1/(1 — t), which 
leads to formula (3), and can be applied to the computation of 7 as discussed above. 
Thb same series (3) can be used to advantage in computing tan~' 1/239 by means 
of the formula 

(4) tan-' ^ - tan-' ^ - 2 tan- ^ . 
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Foreword 

One of the first large electnimc computers ever built, the ENIAC was designed 
and constructed at the Moore fkrhool of Electrical Engineerings University of 
Pennsylvania, Philadelphia^ in 1946, The name stands for "Eleciivnic Numerical 
InUgraior and Calculator"* It was capable of performing 5000 additions per sec^ 
ond and up to 500 multiplications per second* 

Advances and improvements in electronic computers have been unbelievingly 
rapid in the twenty odd years since ENIAC first appeared* Taday*s (1965) machines 
can |ierform 100«000 additions per second and 10,000 muhiplications per second* 
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An ENIAC Determination of tt and e 
to more than 2000 Decimal Places 



Early in June. 1949, Pitifc«K9r JoHN VON NEUMANN expressed an interest in the 
possibility thai the ENIAC might sometime be employed to determine the value 
of TT and e to many decimal places with a view toward obtaining a statistical mea- 
sure uf the randomness of cUstribution of the digits, suggesting the employment of 
one of the formulas: 



Further interest in the project on v was expressed in July by Dr. Nicholas Metro- 
POUS who offered suggestions about programming the calculation. 

Since the possiblity of official time was too remote for consideration, permission 
was obtained to execute these projects during two summer holiday week ends when 
the ENIAC would otherwise stand idle, and the planning and programming of the 
projects was undertaken on an extra-curricular basis by the author. 

The computation of e was completed over the July 4th week end as a practice 
job to gain experience and technique for the more difficult and longer project on 
ir. The reciprocal factorial series was employed: 



The first of the above-mentioned formulas was employed for the computation of 
it; iu advantage over the others will be explained later. The computation of tt was 
completed over the Labor-Day week end through the combined efforts of four mem- 
bers of the ENIAC staff*: Clyde V. Hauff (who checked the programming for tr). 
Miss HoMEf S. McAuJSTER (who checked the programming for eU W. Baskuy 
Fritz and the author, taking turns on eight-hour shifts to keep the ENIAC oper- 
ating continuously throughout the week end. 

While the programming for e is valid for a little over 25(X) decimal places and, 
with minor alterations, can be extended to much greater range, and while the 
programming for v is valid for around 70(X) decimal places, the arbitrarily selected 
Umit of 2000+ was a convenient stopping point for e and about all that could be 
anticipated for a week end's operation for tr. 

While the details of the programming for each project were completely different* 



ir/4 = 4 arctan 1/5 - arctan 1/239 

7r/4 = 8 arctan 1/10 — 4 arctan 1/515 - arctan 1/239 

7r/4 ^ 3 arctan 1/4 + arctan 1/20 + arctan l/I^ 



in conjunction with the GREGORY series 



arctan jc= 2 1)»(2« + l) *^ 
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the grnera! pattern uf pnn edure roughly the same, and both projec ts will be 
dbru^Hed together. In both pmjerts the ENIAC*S divider was employed to deter 
mine a rhoi^en number i of dtgitti of each sureessive term of the series being com- 
puted« the remainder after each division being stored in the ENI A(rS memiiry and 
the digits of eac^h term being added to (or subtracted from) the cumulative total. 
After performing this operation for as many sueressive terms as practicable, the 
remainders for these terms were printed on an I.B.M. card (the standard input- 
output vehicle for the ENIAllK and the pnK*ess was repeated, continuing thn>ugh 
some term beyimd which the digits of and remainders for all further terms would 
be zerns. At this point was printed the cumulative total of the digits of the individual 
terms, which yielded (after adjustment for carry-over) the actual digits <>f the series 
being determined. 

The cards bearing the remainders then were fed into the ENIAd reader, and 
the entire pnK-ess was re|>cated for the next i digits, the ENlAi! reading each 
remainder in turn and placing it before the digits of the appropriate term. Each 
deck of (*ards bearing remainders was then employed to determine the **next*' 
I digits and the **next" deck of **remainder" cards continuing thmugh the firs! 
stfipping point beyond the 2000th decimal place. The cards bearing the cumulative 
totals of sets of i digits of the terms were then adjusted for carry-over into each 
preceding set of i digits. In the case of e this yielded the Bnal resuh; in the case 
of TT ail the above described operations were performed once for each inverse tan- 
gent series, so that ea<'h set of "cumulative total" cards, adjusted for <*:rry-over, 
yielded the digits of ime of the series, the final result being determined by the com- 
bination of these series in appropriate manner. 

The number (»f ptac^es i chosen for each interval of computation, the maximum 
magnitude of each remainder, the amount of memory space available, and the 
details of divider operation (the number of places to which division can be per- 
formed to yield a |N>sitive remainder, and the necessary conditions of relative and 
absolute |Hisitioning of numerator and denominator) ail were interrelated, and 
where oppiirtunity for selection existed, that selection was made which provided 
maximum efficiency of computation. In the case of rr there was imposed the addi- 
tional requirement that identical programming apply for all series employed, and 
for this reason the f<irmula: 

7r/4 = 4 arctan 1/5 - ar< tan 1/239 

was superior to the other two. 

In order to insure abs«tlute digital accuracy, the programming was arranged s( 
that one half applied to c*omputaticm and the other half to checking. Before any 
deck of "remainder'' c ards was employed to determine the next i digits, the cards 
were reversed and employed in the checking sequence to confirm each division 
by a multiplication and each addition by a subtrac tion and vice versa, reproducing 
the previous deck of ••remainder" cards and insuring that the cumulative total 
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reduited to zeru. (In the case of e this was a simple inversion of the computation; 
in the case of ir the factor (2n + 1)^^ in each term made it a more complicated 
affair). After the correctness of each deck was established through this checking, 
the **remainder*' cards were rereversed, and the computation proceeded for the 
next i digits. 

Since the determination of each i digits was not begun until the determination 
uf the previous i digits had been confirmed by checking, the ENIAC stood idle 
during the reversals and rereversaJs and comparisons of the decks in the compu- 
tation of e; in the case of tt, however, the ENIAC was never idle, for operation on 
each series wa» alternated with operation on the other, card«handling on either 
being accomplished while the other was being operated upon by the ENIAC. In 
the case of e, insurance against any undiscovered accidental misalignment of 
cards was provided by rerunning the entire computation without checking, i.e*, 
wilhuul card reversals, confirming the original resuhs; in the case uf ir, the same 
asiiuran(*e was provided by a programmed check upon the identification numbers 
uf each successive card in both computation and checking* 

In the case of e, there wa? printed (in addition to each **remainder" card) a card 
containing the current / digits of (n!)"^ for n « 1, 2, 3 . . . ; in the case 

of 7T only remainder and final total cards were printed. 

The ENIAC determinaliims of both tt and e confirm the K)8— place determina- 
tion of f published in MTAC, v. 2, 1946, p. 69, and the place determination 
of n published in MTAC. v. 2, 1947, p. 245, as corrected in MTAC, v, 3, 194fi, 
p. 1819. 

Only the following minor observation is offered at this time concerning the ran- 
domness of the distribution of the digits. Publication on this subject will, however, 
be forthcoming MMm. A preliminary investigation has indicated that the digits of 
e deviate significantly from randomness (in the sense of staying closer to their 
expectation values than a random sequence of this length normally would) while 
for IT no significant deviations have so far been detected. 

The programming was checked and the first few hundred decimal places of each 
constant were determined on a Sunday before each holiday week end mentioned 
above, the principal effort being made on the longer week end- The actual required 
machine running time for both computation and checking in the case of e was 
around U hours, though card-handling time approximately doubled this, and the 
recomputation without checking added about 6 hours more; actual required ma* 
chine running time (including card^handling time) for tr was around 70 hours. 

The following values of tt and e have been rounded off to 2035D and 2010D 
respectively. 



5r = 3.14159 26535 89793 23846 26433 83279 50288 4I97I 69399 37510 

58209 74944 59230 78164 06286 20899 86280 34825 34211 70679 

82148 08651 32823 06647 09384 46095 50582 23172 53594 08128 

481U 74502 84102 70193 85211 05559 64462 29489 54930 38196 

44288 10975 66593 34461 28475 64825 37867 83165 27120 19091 

45648 56692 34603 48610 45432 66482 13393 60726 02491 41273 

72458 70066 06315 58817 48815 20920 96282 92540 91715 36436 
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78925 90360 01133 05305 48820 

33057 27056 57595 91953 09218 

07446 23799 62749 56735 1885/ 

98336 73362 44065 66430 86021 

60943 70277 05392 17176 29317 

00056 81271 45263 56082 77857 

14684 40901 22495 34301 46549 

42019 95611 21290 21960 86403 

51870 72113 49999 99837 29780 

50244-59455 34690 83026 42522 

71010 00313 78387 52886 58753 

59825 34904 28755 46873 11595 

18577 80532 17122 68066 13001 

38095 25720 10654 85863 27886 

03530 18529 68995 77362 25994 

55748 57242 45415 06959 50829 

81754 63746 49393 19255 06040 

858^6 16035 63707 66010 47101 

94482 55379 774^2 68471 04047 

93313 67702 89891 52104 75216 

25338 24300 35587 64024 74964 

67823 54781 63600 93417 21641 

55706 74983 85054 94588 58692 

32116 53449 87202 75596 02364 

63698 07426 54252 78625 51818 

81647 06001 61452 49192 17321 

16136 11573 52552 13347 57418 

45477 62416 86251 89835 69485 

56887 67179 04946 01653 46680 

82796 79766 81454 10O95 38837 

73929 84896 08412 84886 26945 

06744 27862 20391 94945 04712 

46776 46575 73962 41389 08658 

94657 64078 95126 94683 98552 

2.71828 18284 59045 23536 02874 

95749 66967 62772 40766 30353 

27466 39193 20030 59921 81741 

59563 07381 32328 62794 34907 

15738 34187 93070 21540 89149 

82264 80016 84774 11853 74234 

55170 27618 38606 26133 13845 

67371 13200 70932 87091 27443 

92836 81902 55151 08657 46377 

77078 54499 69967 94686 44549 

77361 78215 42499 92295 76351 

28869 39849 64651 05820 93923 

30123 81970 68416 14059 70198 

53118 02328 78250 98194 55815 

96181 88159 30416 90351 59888 

87922 84998 92086 80582 57492 

84875 60233 62482 70419 78625 

49146 31409 34317 38143 64054 

76839 64243 78140 59271 45635 

01157 47704 17189 86106 87396 

02123 40784 98193 34321 06817 
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46652 13841 46951 94151 16094 

61173 81932 61179 31051 18548 

52724 89122 79381 8*011 94912 

39494 63952 24737 19070 21798 

67523 84674 81846 76694 05132 

71342 75778 96091 73637 17872 

58537 10507 92279 68925 89235 

44181 59813 62977 47713 09960 

49951 05973 17328 16096 31859 

30825 33446 85035 26193 11881 

32083 81420 61717 76691 47303 

62863 88235 37875 93751 95778 

92787 66111 95909 21642 01989 

59361 53381 82796 82303 01952 

13891 24972 17752 83479 13151 

53311 68617 27855 88907 50983 

09277 01671 13900 98488 24012 

81942 95559 61989 46767 83744 

53464 62080 46684 25906 94912 

20569 66024 05803 81501 93511 

73263 91419 92726 04269 92279 

21992 45863 15030 28618 ^9745 

69956 90927 21079 75093 02955 

80665 49911 98818 34797 75356 

41757 46728 90977 77279 38000 

72147 72350 14144 19735 68548 

49468 43852 33239 07394 14333 

56209 92192 22184 27255 02542 

49886 27232 79178 60857 84383 

86360 95068 00642 25125 20511 

60424 19652 85022 21066 11863 

37137 86960 95636 43719 17287 

32645 9958! 33904 78027 59009 

59570 98258 

71352 66249 77572 47093 69995 

54759 45713 82178 52516 64274 

35966 29043 57290 03342 95260 

63233 82988 07531 95251 01901 

93488 41675 09244 76146 06680 

54424 37107 53907 77449 92069 

83000 75204 49338 26560 29760 

74704 72306 96977 20931 01416 

21112 52389 78442 50569 53696 

05987 93163 68892 30098 79312 

48220 82698 95193 66803 51825 

98294 88795 52036 25094 43117 

37679 52068 52823 76464 80429 

30175 67175 61552 06981 12509 

85195 45807 27386 67585 89422 

79610 48419 84445 65463 24496 

20900 21609 90235 30436 99418 

62531 52096 18369 08887 07016 

49061 30310 72085 10583 75051 

96552 12671 54688 95703 50354 

01210 05627 88023 51930 33224 
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74501 58539 04730 41995 77770 93503 66041 69973 29725 08S68 

76966 40355 57071 62268 44716 25607 9^6 51787 13419 51246 

65201 03059 21236 67719 43252 78675 39855 89448 96970 96409 

75459 18569 56380 23637 01621 12047 74272 28364 89613 42251 

64450 78182 44235 29486 36372 14174 02388 93441 24796 35743 

70263 75529 44483 37998 01612 54922 78509 25778 25620 92622 

64832 62779 33386 56648 16277 25164 01910 59004 91644 99828 

93150 56604 72580 27786 31864 15519 56532 44258 69829 46959 

30801 91529 87211 72556 34754 63964 47910 14590 40905 86298 

49679 12874 06870 50489 58586 71747 98546 67757 57320 56812 

88459 20541 33405 39220 00113 78630 09455 60688 16674 00169 

84205 58040 33637 95376 45203 04024 32256 61352 78369 51177 

88386 38744 39662 53224 98506 54995 88623 42818 99707 73327 

61717 83928 03494 65014 34558 89707 19425 86398 77275 47109 

62953 74152 11151 36835 06275 26023 26484 72870 39207 64310 

05958 41166 12054 52970 30236 47254 92966 69381 15137 32275 

36450 98889 03136 02057 24817 65851 18063 03644 28123 14965 

50704 75102 54465 01172 72115 55194 86685 08003 68532 28183 

15219 60037 35625 27944 95158 28418 82947 87610 85263 98139 

55990 06738 

Values of the auxiliary numbers arccot 5 and arcccit 239 to 2035D are in the pos- 
tteHsion uf the author and aim have been deposited in the library of Bmwn Univer- 
sity and the I'MT File' of MTAC. 
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Significantly, the r«al numbera of elementtry algebra fall into two dUjoint sets: 
(1) ihv r<4tiunal numben. and (2) the irrational numbers. A raiionai number in 
a number that can be expretaed as the ratio of two integers, as, for example, 
f • Every rational number when expressed in decimal fracliim form 

yields either a terminating decimal or a repeating decimal. Thus, 

l-«,875, and ^ .181818 ... . 

It is not difficult to show that between any two rational numbers there exist in- 
finitely many other rational numbers. Thus, if the rational numbers were associated 
with p<iint8 on a line, it would seem as if the line were **completely filled*' with 
points. 

Ahhough it is difficult to picture it, such a line is nof completely filled with points. 
Strangely enough, it is full of ''holes", that is, there are many points which have 
no rational numbers assij^ied to them. The numbers that "belong'* to these points 
are called irrational numbers. 

An irrational number is a number that is not rational, that is, it cannot be ex* 
pressed as the quotient of two integers. The existence of non-rational numbers is 
easily shown. A claiisical proof was given by Pythagoras over 2000 years ago, as 
follows. Assume that V¥ is rational. Let vT«= a/6, where a and 6 are relatively 
prime. Then 

^ \bl b^ * 
or 263«a^ (1) 

Hence a* is an even number; therefore a is also an even number. But if a is an even 
number* it can be expressed as 2k. where * is any positive integer. Thus 

or l^^W. 

Hence 6* is an even number, and therefore b is an even number. Since both a and 
6 have been proved to be even numbers, the assumption that a and b are relatively 
prime is false, and so the assumption that VT »= a/6 is false. In short, cannot 
be expressed as the ratio of two integers. The prot>f can be generalized to \//V , 
where N is any integer which is not the square of another integer. 

One of the reasons for the many attempts to find the value of ir to so many deci- 
mal places is the desire to learn something about the distribution of the digits in 
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the extended approxinuition of It b«» been proved th«t «■ is an irrationaJ number, 
that ii a number which when expreiiied as a decimal in base 10. yields a non- 
terminatinfi. nun-repeating decimal. 

An irrational number is said to be a normal number if all the digiu occur with 
equal frequency, and if all blocks of digits of the same length occur with equal 
frequency, from the standpoint of the theory of numbers and higher analysis, 
mathematicians are curious about the distribution of the digits in the numerical 
approximation of ir. It is believed that ir is a normal number with respect to base 
10, but it is not yet known whether v is normal to any base. These and related 
questions are of consideralJ.' interest to modern mathematicians. 
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The Evolution of Extended Decimal 
Approximations to tr 

by J. W. Wkench« Jit, 

In hb hUtoHcai survey uf the cUmic problem of ^'squaring the circle/* Profet- 
•or E. W, Hobson [l]*difttiiiguished three distinct periods, characterized by funda- 
mental diiferetices in method, immediate aims, and available mathematical tools. 

The first period— the so-caUed geometrical period --extended fnim the earliest 
empirical determinations of tli^ ratio of the circumference of a circle to its diameter 
to the invention of the calculus about the middle of the seventeenth century. The 
main effort was directed toward the approximation of this ratio by the calculation 
of perimeters or areas oi regular inscribed and circumscribed polygons. 

The second period began in the middle of the seventeenth century and lasted 
for more than a hundred yean. During this period the methods of the calculus 
were employed in the development of analytical expressions for ^ in the form of 
infinite series, products, and continued fractions. 

The third period, which extended from the middle of the eighteenth century to 
nearly the end of the nineteenth century, was devoted to studies of the nature of 
the number rr. J. H. Lambert [2] proved the irrationality of ir in 1761, and F. 
Lindemann [3] first established its transcendence in 18S2. 

This article is concerned with the second period and its sequel, which extends 
to the present day. 

According to Hobson [1], the first analytical expression discovered in th^s 
period is the infinite product 

fr^2 2 4 4 6 6 8 8 
2*T'3'3'5'5*7'7'9"' * 

which was published by John WalUs [4] in 1655. 

Lord Brouncker, the first president of the Royal Society, about 1658 found the 
infinite continued fraction 

^-J- JLiLil 

4"l + 2+2+2+'" • 
which was shown subsequently by Euler to be equivalent to the alternating series 

4 ^ l^S^7^9 
known to C. W. Leibnia in 1674 

^NMbm im bmekcu mlef to Om rwfeftmom iUlMi mt end iif the ftitkle. 
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The great nu^jority of calcalatioiu of ir to m*ny decimal placet have been based 
upon the power »erie> 

orcton x — x — ^ + 'g-— ... ,— 1 S xS 1, 

which was discovered in 1671 by jonies Gregory [5]. He failed, however, lo note 
explicitly the special case corresponding to x >- 1. which is ascribed to Leibniz. 
Sir Isaac Newton [6] in 1676 discovered the power series 

^lx»^l'3x»^ 
arc.mx-=x + r- + 5:j-^.... 



- 1 S X S 1, 

which has been used by a few computers of n. 
In 1755 Leonhard Euler [7] obtained the following useful series: 

-rctanx»y^{l+|(y^) + |^(y^)V...}. 

It was by means ofGregory's series, taking x - 1/VF, that Abraham Sharp [8], 
at the suggestion of the English astronomer Edmund H alley, computed v to 72 
decimal places in 1699, thereby nearly doubling the ^«aicst accuracy (39 ctecimal 
places) attained by earlier computers, who had used geometrical methods. Sharp's 
calculation was extended by Fautet de Lagny [9] in 1719 to 127 decimals (the 1 13th 
place has a unit error). 

Newton set X " — i in his series, and thereby computed «r to 14 places. A Japanese 
computer. Matsunaga Ryohitsu [10], used the same procedure to evaluate ir cor- 
rect to 49 (^imal pliu^s in 1739. About 1^ a Chinese. Chu Hung, calculated 
tr to 40 places (25 correct) by this series [10]. 

Most computers of fr in modem tin^s have used Gregory's series in coiyunction 
with certain arctangent relations. Only nine of these relations have been employed 
to any extent in such computations. We shall now oonsiiter these formulas, ar- 
ranged according to the increasing precision of the approximations computed by 
their use. 

ir 1 3 

1.7-5 arctan r + 2 arctan =r 

4 7 79 

EuIer [7] in 1755 used this relation in ooi\junction \ *h his series for arctan x to 
compute v correct to 20 decimal plwes in one hour. B ^n Georg von Vega [ U ] in 
1794 employed Gregory's series and the preceding rciution to evaluate tr to 140 
decimal places, of which the first 136 were correct. This precision was exceeded 
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1^ that atuiried by an unknowo cakulaior whose manuscript, containing an up- 
imsiximation correct to 1S2 pia^^, was sMn in the Riuicliffe Libraiy at Oxford 
toward the close the eighteenth century, 

II, 7*4 arctan \ — arctan ^ arctan ~ 

Euler published this reiatibn in 1764. It was used by William Rutherford [12] 
in I84f 10 compute to 208 places (152 correct). 

IT 1 1 1 

III. 7 arctan -z + arctan r + arctan r 

4 ^ b a 

This formula was supplied the calculating prodigy Zacharias Dahse [13] by 
Schulz von Strassnitzky of Vienna, Within a period of two months in 1844« 
Dahse thereby evaluated ir correct to 200 places. 

IV, ^ « arctan | + arctan | 

First published by Clmrles Hutton [14] in 1776, this relation was used by W. Leh- 
mann [15] of Potsdam to compute ^ to 261 decimals in 1853. Tseng Chi-hung [16] 
in 1S77 used the same formula to evaluate ^ to 100 decimals in a little more than 
a month. 

7T 11 

V. 7 == 2 arctan r + arctan 



The relation was also pubiiihed by Hutton [14] in 1776, and independently by 
Euier in 1779. Vega [17] u^ it in 1789 to compute 143 decimals (126 crnrect). In 
order to remove the ur^ertainty causoi by ihit discr<:pant approximations of 
Rutherford and Dahse. Thomas Clausen [18] extended the calculation to 248 cor- 
rect decimals in 1847, and Lehmann [15] reached 261 decimals in 1853 by this 
formula, confirming his inctependent calculation of it to the same extent by rela* 
tion Edgar Frisby [19] in Wishington, D. C. used illation V in conjunction 
with Euler*s series to compute ^ to 30 places in 1872. 

ir 1 1 1 

VI. - « 3 arctan ^ + arctan + arctan 

This formula was published by S. L. Loney [20] in 1893, by Carl Stonror [21] 
in 1896. and was rediscovered by R. W. Morris [22] in 1944. By means of this 
formula D. R Feiguson, then of the Rc^al Naval Colle^, Eaton, Chester, England, 
performed a bnghand calculation of tr to 530 decimal places between May 1944 
and May 1945* At that time he discovered a discrepancy between his approxima- 
tion and t}i^ final result of William Shanks -discussed under formula DC -begin- 
ning with the 528th place. The first notice of an error in Shank's well-known 
apfHTOumation appeared in a note [22] published by Feiguson in March 1946. He 
continue his calculatkin of n and in July 1946 published [23] a oorrectk)n to 
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Shank*! value through the 620th deciniai pl^e. Subsequently, Ferguson used a 
dwk calculaiur to reach 710 decimals [24] by January 1947. and finally 808 deci- 
mailt [25] by September 1947. 

VII. |»8arctan~-arctan~-4arctang|^ 

S. Klingenstierna discovered this relatkin in 1730; it was rediscovered more 
than a century later by Schellbach [26]. It was used by C. C. Camp [27] in 1926 to 
evaluate ?r/4 to 56 places. D. H. Lehmer [28] reoomn^nded it in coiuunctwn with 
the next formula for the calculatwn of « to many figures. G. E. Felton on March 3 1 . 
1957 completed a calculatnn of ir to 10Q21 places on a Pegasus computer at the 
Ferranti Computer Cenure in London. This required 33 hours of computer time. 
The result was published to 10000 places [29]. A check calculatwn using formula 
VIII revealed that, because of a machine error, this result was incorrect after 7480 
(tecimal places. 

Gauss [30] investigated the derivation of arctangent relations and reduced it 
to a problem in Diophantine analysis. Relation VIII is one of several formulas he 
developed. J. P. BaUantine [31] substantiated Lehmer's claim that this formula 
is especially effective for extensive calculation, by discussins its use in conjunc- 
tion with Euler's series for the arctangent. 

VIII. I = 12 arctan ^ + 8 arctan ^ - 5 arctan ^ 

Felton carried out a second calculation to 10021 places, and by March 1, 1958 had 
removed all discrepancies from his results, so that the approximations computed 
from formulas VII and VIII agreed to within 3 units in the 10021st decimal place. 
The corrected result remains unpublished. 

IV 4 I 1 

lA. - = 4 arctan ^ - arctan — 

This is the most celebrated of ail the relations of this kind. John Machin, its 
discoverer, computed n correct to 100 decimals by means of it in conjunction with 
{;regory'8 series, and the result [32] appeared in 1706. Clausen [18] in 1847 used 
this relation in addition to Hutton's formula V to compute w to 248 decimal places, 
as has already been noted. 

Rutherford resumed his calculation of v in 1852. using Machin's formula this 
time, as did his former pupil William Shanks. Shanks's first publinhed approxi- 
mation to n contained SSO decimal places, and was incorporated in Rutherford's 
note [33]. published in 1853, which set forth his approximation to 441 decimals. 
Later that year Shanks published his book [34] containing an approximation to 
607 places and isiving all details of the calculation to 530 places. It is now known 
that Shanks's value was incorrectly calculated beyond 527 decimal places. The 
accuracy nf that value was further vitiated by a blunder committed by Shanks in 
correcting his copy prior to publication, with the resuh that similar errors appear 
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in decimal piacen 460^462 and 513—515. These errorti perni&t in Shanks't^ finit 
paper of 1873 [35] containing the exteniion to 707 decimals of his earlier approxi- 
mation. Hi» ici^und paper of thai year [36] which contained his final approximation 
to iTt ipve;^ corrections of theae errors; however, there appears an inadvertent 
typographical error in the 326th decimal place of his final value* In retrospect, we 
now realize that Shanks's first value published in 1853 was the must accurate he 
ever published. 

The accuracy of Shanks^s approximation to at least 500 decimals was confirmed 
by the independent calculations of Professor Richter [37] of Elbing« Germany, 
who in 1853—1854 computed successive approximations to 330, 400, and 500 
places. Richter*s communications do not reveal the formula that he used. 

Machines formula was used by H. S. Uhler in an unpublished computation cor- 
rect to 282 places, which was completed in August 1900. 

F. J. Duarte computed n correct to 200 places by this method in 1902. The re- 
sult was published [38] six years later. 

As a by-product of his calculation of the natural logarithms of small primes, Uhier 
in 1940 noted [39] confirmation to 333 decimal places of Shanks's approximation. 

In December 1945. Professor R. C Archibald suggested that the writer under* 
take the coinputatiim of n by Machines formula in order to provide an independent 
check of the accuracy of Ferguson*s calculations. With the collaboration of Levi 
B* Smith, who evaluated arctan 1/239 to 820 decimal places, the writer computed 
IT to 818 places by February 1947, using a desk calculator. The result was published 
[24] to 808 places in April 1947. and was verified to 710 places by Ferguson in 
a note published concurrently [24]. The limit of 80S decimals in the published 
value was chosen to provide precision comparable to that obtained by P. Pedersen 
[40] in his approximaticm to e. 

Collation Aii this 808-pIace approximation with results obtained by Ferguson 
later that year revealed several ern>neous figures beyond the 723rd place in the 
writer's approximatiim to arctan Vs. These errors vitiated the corresponding figures 
in the approximation to tt. Corrections of these errora and extensions of Ferguson's 
results appeared in a joint paper [25] by Ferguson and the writer in January 1948, 
which concluded with an 808*piace approximation to rr of guaranteed accuracy. 

Subsequently, Smith and the writer resumed their calculations and by June 1949 
had obtained an approximation to about 1120 decimal places [41]. Before final 
checking of this extension could be completed, the ENIAC (Electronic Numerical 
Integrator and Computer) at the Ballistic Research Laboratories* Aberdeen Prov- 
ing Ground, was empkiyed by George W. Reitwte^ner and his associates in Sep- 
tember 1949 to evaluate rr to about 2037 places (2040 working decimals) in a total 
time (including card handling) of 70 hours [42]. Machines formula was also used 
in this computation. 

In November 1954. Smith and the writer extended their calculation to 1150 
places, and in January 1956 reverted to this work once more to attain their final 
result, which was terminated at 1160 places, of which the first 1157 agree with 
those obtained on the ENIAC. 

2S 



ERIC 



20 



A calculation of n was performed in duplicate on the NORC (Naval Ordnance 
Research Calculator) tn November 1954 and in January 1955 aa a demonstration 
problem, prior to the delivery of that computer to the U. S. Naval Proving Grounds 
at Dahlgren, Virginia. Again, Machines formula was selected, and the calculation 
was completed to 3093 decimal places in 13 mihutes running time, A report of this 
work, in which the value of ^ was presented unrounded to 3089 decimal places, 
was published by C. Nicholson and J. Jeenel [43] of the Watson Scientific Com- 
puting Laboratory, in New York. 

In January 1958, Francois Genuys [44] programmed and carried out the evalua- 
tion of TT correct to 10000 decimal places on an IBM 704 Electronic Data Process* 
ing System at the Paris Data Processing System at the Paris Data Processing Cen- 
ter. Machines formula in conjunction with Gregory's series was used. Only 40 
seconds were required to attain the 707 decimal-place precision reached by 
Shanks, and one hour and forty minutes was required to reach the lOOOO places 
of the final result. 

On July 20, 1959, the program of Genuys was used on an IBM 704 system at the 
Commissariat a TEncrgie Atomique in Paris to compute ^ to 16167 decimal places. 
This latest approximation Is unpublished at present. 

TABLE 1 

Cumulative distribution of the first 16000 decimal digits of tt 



Digit 

Thousand - 





0 


I 


2 


3 


4 


5 


6 


7 


8 


9 


I 


93 


116 


103 


102 


93 


97 


94 


95 


101 


106 


2 


1 82 


212 


207 


ISS 


195 


205 


200 


197 


202 


212 


3 


259 


309 


303 


265 


318 


315 


302 


287 


310 


332 


4 


362 


429 


408 


568 


405 


417 


398 


377 


405 


431 


5 


466 


532 


496 


459 


508 


525 


513 


488 


492 


512 


6 


557 


626 


594 


572 


613 


622 


619 


606 


582 


609 


7 


657 


733 


692 


686 


702 


730 


708 


694 


680 


718 


8 


754 


833 


SU 


781 


809 


834 


816 


786 


764 


812 


9 


855 


936 


911 


884 


910 


933 


914 


883 


854 


920 


10 


968 


1026 


1021 


974 


1012 


1046 


1021 


970 


948 


1014 


n 


1070 


1099 


nil 


1080 


1133 


1150 


1129 


1070 


1031 


1127 


12 


1162 


1193 


1214 


1176 


1233 


1262 


1227 


1166 


1144 


1223 


13 


1266 


1314 


1316 


1272 


1343 


1358 


1324 


1260 


1246 


1301 


14 


1365 


1416 


1419 


1383 


1440 


1455 


1426 


1344 


1339 


1413 


15 


1456 


1513 


1511 


1491 


1553 


1549 


1520 


1441 


1458 


1508 


16 


1556 


1601 


1593 


1602 


1670 


1659 


1615 


1548 


1546 


1610 



The motivatiun of modern calculations of n to many decimal places was conjec- 
tured by Professor P. S. Jones [45] in 1950 as being attributable to "intellectual 
curiosity and the challenge of an unchecked and long untouched computation." 
This reason for undertaking such work should be supplemented by reference to 
the recurrent interest in determining a statistical measure of the randomness of 
distribution of the digits in the decimal representation of ir. 
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Augusta* Dc Morgan [46] drew attention to the deficiency in ihc number of 
appearance* of the digit 7 in Shank»*« «)7.place approximation to In 1897 E. B, 
Escott [47] raised the quculion whether the de6ciency of 7*8 noted in Shanks's 
6nal approximation could be explained. 

In June 194S, the late- Professor John von Neumann expressed an interest in 
utilizing the ENIAC to determine the value of tt and e to many places as the basiti 
for a statiftitral study of the distribution of their decimal digits. A statistical treat- 
ment of the first 2000 decimal digits of both rr and e was published by N* C, Metro- 
polis, Reitwiesner, and J. von Neumann [48], Furtheranalysisof these data was 
performed by R. E, Greenwood [49], using the coupon coUector*s test. A count of 
each of the decimal digits appearing in the NORC approximation appears in the 
paper of Nicholson and Jeenel [43]. A number of recent investigators have dis- 
cussed the distribution of digits in Shanks*s approximation and in the corrected 
value of IT. These investigators include F. Bukovszky [50], W. Hope-Jones [51], 
E. H. Neville [52], and B. C. Brookes [53], 

The writer has recently completed a count by centuries of the 16167 decimal 
digits constituting the fractional part of the latest approximation to ir. An abridg- 
ment of this information is presented in the accompanying table. 

The standard }^ test for goodness of fit reveals no abnormal behavior in the dis- 
tribution of digits in this sample; in particular, there appears to be no basis for 
supposing that tr is not simply normal [54] in the decimal scale of notation. It has 
been pointed out recently by Ivan Niven [55] that the normality of such numbers 
as e, and VT^ has yet to be proved. 

Numerical studies directed toward the empirical investigation of the normality 
Iff TT clearly require inc reasingly higher decimal approximations, which can best 
be obtained by use of ultra-high speed electronic computers now under design and 
development. 
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Did You Know That w Has Been Calculated 
to 100^65 Decimal Places? 

Daitiel ShAfiki and John W. Wrench, Jr., both of the David Taylor Model Basin 
in Waihington, D.Cm have calculated the values of n and e to lOOJ^D on an IBM 
7090 iyatenu' The computation of v was performed July 29, 1961 at the IBM 
Datacenter in New York and required 8 hours 43 minutes while the evaluation of 
e required 2.5 hours. 

Before going into the question of M;Ay such calculations are made, let*s backtrack 
a bit and see what has been done in the past. 

The Bible is content with a value of 3 for the ratio of the circumference of a circle 
to its diameter but Archimedes was able to assign limits to tr between 3V? to 3^^/7i. 
The Egyptians had managed to evaluate, somehow, ir as about 3.16 while the Baby- 
lonians used the same value, 3, as the Bible. 

It is remarkable that the Chinese Astronomer Tsu Ch*ung-Chih discovered a 
simple fraction in the 5th Century that gives the value of ir accurate to six decimal 
places: 




3.1415929 . . . 



From the middle of the 17th century many approximation expressions in the form 
of infinite series of one kind or another were developed. Evaluations of tr to as many 
decimal places as the patience of the computer could stand followed rapidly* 

tt was computed to 72D by Abraham Sharp in 1699; to 127D by Fautet de Lagny 
in 1719; to 49D by the Japanese computer Matsunaga Ryohitsu using a {X)wer 
series developed by Sir Isaac Newton; to 140D in 1794 by Baron Georg von Vega 
(but only his first 136D were correct); to 40D by Chu Hung in 1800; to 152D by an 
unknown computer at the close of the 18th centiu^; to 208D by William Ruther- 
ford in 1841, using one of Euler^s arctangent relations; to 261 D (twice by different 
methods) by W. Lehmann in 1853.* 

The most celebrated calculation of tr was made to 7p7D by William Shanks on 
and off for the 20*year period from 1853 to 1873* It was not until 1945 that Shanks 
was found, by D. F* Ferguson, to have erred at the 528th decimal place. 

J. W* Wrench, Jr. and D* F. Ferguson calculated n to 808D in 1947 to match the 
evaluation of e at that time. 

th0««ii t» w ciUcsUmi c#ji tw gftinrd hj rcmsidirriaf ttuif. in O^mMuy* ihe vaJi^ of « to 35 dr^iniAi 
pUcm k c^M tkf UdojphiM nuaiber in memofj of t.^4o{>h vas CcuWp, • Grrmsa luAtbrmttirUn. Vaq Ceulrn, in 
IMy cdc^dMd w te SSD Md rcqtmlvd ihaM t^ft vntar be le«cHM on bit tomb^onr «• ah epiUph. He died At Ibr 
^ «r 79 Awi ffiwibAlnnf was dutifuUf ktftcf^cd As requetled. 
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All thr above calculatbnn were dune longhand (including Shanks 707DU or with 
a deiik calculator. Subsequently, electronic computeni were used and extended 
IT evaluations followed: to 1I20D in June 1949; to 2037D in September 1949 (taking 
70 hourti); to 30930 in November 1954 and January 19S5 (taking only 13 minutes): 
to 1(X}60D (in 1 hour 40 minutes) in January 1938 by Francois Genuys on an IBM 
704 Ellectronjc Data Processing System in Paris; and, almost finally, to 16167D in 
July 1959. 

The latest calculation is that mentioned in the first paragraph. 

Simon Newcomb, the astronomer and mathematician, once remarked about tr 
that ten decimal places would suffice to give the circumference of the earth ac- 
curate to a fraction of an inch and that thirty decimal places would give the circum- 
ference of the known universe to microscopic accuracy! 

Why in the world is such apparently pointless work being dune? 

One practical reason is that new computers can be checked by programming 
problems with known answers. 

The more interesting reason — more interesting to recreational mathematicians, 
anyway — is to tind out, by actual calculation, whether such numbers as tt, e or 
VT are "normal" numbers. That is, whether the digits 1, 2, 3, 4, 5, 6, 7, 8, 9, 0 
^cur in a statistically random distribution — do they each cn-cur approximately 
10% of the time? 

A count of the first 1600DD of tt reveals no abnormal distribution, n appears to 
be a "normal" number.^ At present there is no proof of the normality of such num- 
bers. It is not even known if the ten consecutive digits 1234567890 occur at least 
once in the infinite decimal evaluation of ir. 

Shanks and Wrench estimate that computers will become available, in the next 
5 to 7 years, which will be able to calculate n to 1,000,000 decimal places. (The 
IBM 7090 which performed the feat to 100,265D in 8 hours 43 minutes would re* 
quire months to do the calculation to I,000,OOOD.) 
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Epilogue 

We have ubiirrved that in general, two methods of thinking have bt^en employed 
in vumputing the value of tt: (1) the geometric approach* and (2) the anolytieal 
approaeh* It should be mtted that the man who used the first of these methods 
thought of TT OS equivalent lu a geometrieal ratio, even as the (^reek geometers 
considered the ratio of two line segments when studying metric properties of 
geometric figures. On the other hand, mathematicians using the second method 
think of ?r nut as the ratio of two lengths, but as the symbol for a specific number 
(like the number e 2.718 . . . ) which enters many fields of mathematical analysis 
from theoretical considerations rather than from any question of practical measure* 
ment. In this connection, one of the most remarkable of mathematical relations is 
that which associates rr and e, namely, e*^ + 1 ^ 0. The number e is itself a unique 
constant, being the limit of the expansion (1 + l/n)" as n increases without limit. 
We know that e is not «mly an irrational number, but, like tt, it is also a transcen- 
dental number, that is. a number which is not the riMit of a polymimial equatitm with 
ration al co efficients. The number i is the pure imaginary unit, where =-^1, or 
i = V-- 1. That the product of { and ir, applied to e as an exfmnent, should yield 
the simple integer ^ 1, is indeed an amazing relation. 
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